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I. INTRODUCTION 



The charged current interaction between neutrinos and their associated leptons in medium give an effective matter 
potential which can lead to a resonant flavour conversion called the Mikhcyev-Smirnov- Wolfcnstcin (MSW) effect [Hi]. 
Such behaviour solves elegantly the solar neutrino deficit problem first pointed out by pioneering experiment of R. 
Davis in the 60's Such phenomenon is of crucial importance for the propagation of neutrinos in the supernova 
environment whose flux detection could yield precious information concerning the dynamics of the density profile or 
fundamental neutrinos properties like the hierarchy or the value of the third mixing angle. In supernovae, depending 
on the mass hierarchy, the electron neutrino may encounter one or two resonances via charged current, while the 
muon or tau neutrinos will only interact via neutral current, indistinguishably at the tree level in the Standard Model 
(SM). Indeed, due to the absence of muon or tau particle in such environment, coherent forward scattering may only 
intervene via neutral current to which all flavours are sensitive. Considering that matter interaction can also be seen 
as neutrino index of refraction [1, HI, Botella et al. [||, after showing that correction at O (a) were negligible in the 
case of a neutral medium Q, proved that differences at one- loop in the neutrino index of refraction could arise for 

muon and tau neutrinos at order O ( nBi "i g w Later, supersymmetric radiative corrections, in the case of the 

minimal supersymmetric model (MSSM) has been partially calculated [|[ showing that it could give potentially much 
larger radiative effects than in the SM. It is interesting to note that such radiative corrections have been for a longtime 
considered as negligible or without observable consequences for the supernova environment. Such consideration was 
probably true when not taking into account the neutrino-neutrino interaction which has dramatically changed the 
vision we have of neutrino pro pagatin g in an exploding star and which h as g one through an intense investigation 
S M OS HI 01 Q EE Htl M MM M H M M M M M M M MM M M ■ Actually few papers have 
shown the importance of the one-loop correction in addition with the neutrino-neutrino interaction. First it has been 
shown that for early time after post-bounce when the matter density profile is very high, neutrinos can encounter 
the n — t resonance and possibly modify the v e flux [25L H3- Unfortunately, it has also been shown that in such 
case the high density in addition with a multi-angle neutrino-neutrino interaction would make those effects vanish 
[3lj . Besides those works, only one paper 1 [2^ has shown the importance of V^ T via the influence of the CP-violating 
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1 In [27j . it was noted that the CP- violating phase could influence sensibly the fi and r neutrino flavor. 
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phase 5 contained in the MNSP matrix and whose value is still unknown. Such term induces effects of a non-zero 
CP-violating phase on the electron neutrino fluxes inside and outside the supernova. 

The goal of this paper is to calculate such corrections in the SUSY framework with and without taking into account 
R-parity breaking interactions. It is organized as follows: Sec. 2 introduces the theoretical framework where we 
briefly remember how to calculate radiative corrections in this context and more importantly where we introduce the 
supersymmetric framework. The corrections with R-parity conservation are calculated in Sec. 3 and Sec. 4 is dedicated 
to the case where the R-parity is broken. Before concluding, we explicitly demonstrate in Sec. 5 the influence of V^ T 
on the v e flux when the CP-violating phase is non-zero. 



II. THEORETICAL FRAMEWORK 



A. The calculations of the radiative corrections 



Neutrinos interaction through matter can be described using indices of refraction and in this case, the evolution 
equation of neutrinos with matter, omitting the neutrino-neutrino interaction is: 




1 



-U 



'"'12 u 






Am| 2 



Vv 



An efl 















1 y 


(S) 








An TM J _ 





(1) 



where U is the MNSP matrix, p v the neutrino momentum, Am| = m 2 . — rn~. and An a p = n Ua — n U/3 . 

To size the effect of such radiative correction on the neutrino propagation, we study the one-loop effect on the 
scattering amplitude matrix which describes the interaction of neutrinos with matter: 



(2) 



In the calculations we will use the same approximations as in |6| and [8(] i.e, neutrinos are propagating though an 
unpolarized medium at rest. Consequently, the neutrino index of refraction can be written as: 



Pv(n Vl - 1) 



-V2G F 

f=u,d, 



(3) 



where Nf is the number density of fermion / in the medium. Therefore, the interesting parameter to study in order 
to size the radiative correction to matter interaction is which is defined at the tree-level by 



rtv 



T 3 (f L )-2Q f sin 2 6 w + 5, 



(4) 



with sin" 8 



w 



/m 2 z ~ 0.23, while Qf and T^(fi,) are respectively the electric charge and the third component 



of the weak isospin of the fermion fc. To parametrize the loop corrections to Botella et al. 
some kind of an arbitrary but convenient way a new by 



have defined in 



r 3 (/ L )-2Q / A^s 2 - 
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(5) 



for / 7^ I where the p V1 ^ includes the f-dependent (box) diagrams contributions. Since the \ ut J are chosen to be 
independent of the /, in a electrically neutral medium, they will not contribute to An Tfi . Consequently, An TfJi will 
only be sensitive to Apt = p VT ' — p v ^$ . Note that we are only interested about the difference between n v and n Vr 
indices of refraction because the loop correction to An ell = will be negligible since electron neutrinos 

already encounter charged current interactions with matter. In the SM, the correction have been calculated and are 
found to be small : 



An T ^ = V„ T = eV e ~ 5.4 x 10~ 5 V e (6) 

where V^ T is the effective matter potential to tau neutrinos due to one-loop corrections and e the ratio between the 
and V e which yields the size of the loop correction in comparison with the charged-current matter potential for 
electron neutrinos V e . 
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B. The supersymmetric framework 



The most significant theoretical issues of the Standard Model (SM) are the hierarchy problem for the Higgs mass and 
the non-unification of the gauge couplings . Supcrsymmctry allows us to address these problems and is an attractive 
candidate for new physics beyond the Standard Model (SM). In this theory a supersymmetric particle called LSP 
(Lightest Supersymmetric Particle) is a natural candidate for dark matter. Among various supersymmetric models 
the most extensively studied is the minimal supersymmetric model (MSSM). 

The MSSM contains the minimum number of fields to describe the known SM particles and their superpartners. 
These fields can be gathered into chiral superfiels and vector superfields. A chiral superfield $ is a multiplet which 
contains a scalar field (z), a fermionic field (ifi) and an auxiliary field (F): $ = (z, ip, F). A vector superfield V 
is a multiplet which contains a bosonic field (v^), a fermionic field (A) and an auxiliary field (D): V = (v^, A, D). 
F and D are auxiliary fields, they do not have kinetic terms. They are eliminated by the minimization equations 
of the Lagragian. To SM fermionic fields (leptons, quarks), we have supersymmetric scalar fields (sleptons, squarks) 
associated. Concerning SM vector bosons (t/(l)y, SU(2)l, SU(3)c gauge bosons), fermionic fields called gauginos 
(bino, winos, gluinos) are associated to. Finally we have fermionic fields called higgsinos associated to Higgs bosons. 
Higgsinos and gauginos will mix to generate ncutralinos and charginos. 

A superfield is a function of spacetime coordinates and so-called superspace coordinates which appear as anticom- 
muting Grassman variables (34|. It can be expanded in terms of its component fields and these Grassman variables. 
Products of superfields can be developed resulting in products of individual particle and sparticle fields. 

We call supcrpotcntial all the rcnormalizable products of chiral superfields in a given supersymmetric model. The 
MSSM's supcrpotcntial is: 



W mssm — h t Q.H u T R — h b Q.H d B c R — h T L.H d L c R 
+ fiH u .H d 



(7) 



where H U ,Q... are chiral superfields. The fermionic part of the supersymmetric lagrangian is obtained through 
a procedure where the superfields of the superpotential are expanded in terms of component fields and then the 
superspace coordinates of the result are integrated out: 

CmSSM = h^Q-HutpT^-hblpQ.HdtpB^-hrlpL-HdlpL^ ,o\ 

+H^h u 4h 4 + ...(+h.c) 1 ' 

In Supcrsymmctry, we need two Higgs bosons to give masses to the other particles. 
Despite its simplicity, there are two unexplained hierarchies, within the MSSM: 

- the so-called //-problem [35[. It arises from the presence of a mass ^-term for the Higgs fields in the super- 
potential. The only two theoretical natural values for this parameter are either zero or the Planck energy 
scale. However, we need fi > lQQGeV to satisfy LEP constraints on the chargino masses and fi < MguSY 
for a destabilization of the Higgs potential in order to have non-vanishing v.c.v. for the scalar Higgs fields 
(< H u , d >^ 0) 

- The other hierarchy with an unknown origin is the one existing between the small neutrino masses (smaller than 
the eV scale) and the electroweak symmetry breaking scale (~ lOOGeT^). 

In this paper, our framework is the Next-to Minimal Supersymmetric extension of the Standard Model (NMSSM) 
[3(| . This model provides an elegant solution to the /i-problem via the introduction of a new gauge-singlet superfield 
S that acquires naturally a v.e.v. x of the order of the supersymmetry breaking scale, generating an effective fx 
parameter (Xx — Li e g) of order of the electroweak scale. Furthermore, this model explains the second hierarchy by 
generating two neutrino masses at tree level through R-parity breaking [37l | as we will see below. 

The NMSSM's superpotential is: 

Wnmssm — h t Q.H u T R — hbQ.HdBft — h T L.HdL R , , 

+ XSH u .Hd + § S 3 . { ' 

Then, the fermionic part of the NMSSM lagrangian is: 

CnMSSM = ht^Q.HutpT^-hblpQ.HdipB^-hrlpL-HdlpL^ / 1fl -, 

+\S*l> Hu 4 Hd + nS^sips + ■■■ i+h.c) [W > 



The NMSSM contains the following particles: 
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. Standard Model fcrmions and their left and right scalar supersymmctric partners (sfcrmions) 

. Standard Model gauge bosons 

. 3 neutral scalar Higgs (hi, hi-, h^) 

. 2 neutral pseudo-scalar Higgs (ai, 02) 

. lcharged Higgs (H ± ) 

. 2 charginos (X12) originally from mixing between charged fermionic superpartners of the gauge bosons (gauginos) 
and charged fermionic superpartners of Higgs bosons (higgsinos) 

. 5 neutralinos (xl 5) of which the lighest called lighest supersymmctric particle (LSP) is gererally stable and is 
thence a natural candidate for Dark Matter. They come from the mixing between neutral gauginos et neutral 
higgsinos. 

The NMSSM phenomenology contains the MSSM phenomenology if we take A — > 0, k — > by keeping /x e ff of the 
order of My/. 

Experimental limits on supersymmctric particle masses obviously show that SUSY has to be broken. Thus, the 
supersymmctric particle masses will be different from Standard Model particle masses. The soft SUSY breaking terms 
in the NMSSM are: 

a) mass terms for scalar particles: 

m\ u \H U \ 2 + m\ d \H d \ 2 + m 2 s \S\ 2 + m 2 Q \Q z \ 2 + m 2 jT R \ 2 + m 2 jB R \ 2 + m\ 3 |L 3 | 2 + m 2 R \L R \ 2 + ... 
where all the fields are scalar fields. 

b) mass terms for gauginos: 



c) soft terms associated to the superpotential: 



XA x SH u .H d + -A K S 3 + h t A t Q.H u T R - h b A b Q.H d B c R - h T A T L.H d L c R + ... (+h.c) 



where all fields arc scalar fields. 



III. R- PARITY CONSERVED SUSY CORRECTIONS 



In Supersymmetry, one usually considers all sfermions to be exactly degenerate at the GUT scale, and obtain their 
low energy splittings from the renormalization group evolution of the soft parameters and from terms arising after the 
electrowcak symmetry breaking. In this way, although squarks become significantly split from slcptons, the splittings 
among the masses of different slcpton generations arc only due to the small r- Yukawa coupling. This usually implies 
that m? — to? is 0(m 2 .), and hence the radiative effects on the v U T indices of refraction are, in this case, not larger 
than the SM ones. In the following, we shall consider in the following a large tl~-t r mixing. This occurs for large 
values of the effective [i parameter and tan/3, or for large values of the parameter A of the trilinear soft terms, in 
which case the splitting can be 0[m T (A T + fi c gtanf3)] as we can see below. 

We give here the stau and smuon mass-squared matrices where we neglect splittings due to D-terms among charged 
and neutral sleptons or among £l and £ R . The diagonal terms are respectively the £ R and £l mass-squared terms. 

For tau slcptons: 

_( m 2 T + m\ R m T (A T + /i eff tan/?) 
r I m T (A T + ^ e fftan/3) m 2 . + m\ 3 

and for muon sleptons: 





where 
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~ , m 2 T and m? are of the order of M?„ 



I' n 



susy 



0.1 - ITeV 2 , 



(A T + n c f[tan(3), (A u + (j, e f[tanj3) are of the order of M, 
m T = 1.8GeV, m u = 105MeV. 



susy : 



In the 2 nd generation case, we can neglect the mixing between fin and fiL because the off-diagonal terms are 
negligible w.r.t. the diagonal terms. But in the third generation, we have to consider fi, T2 instead of tr, t£. 
The SUSY contribution to An Tfl can then be larger than the SM one. 

In the following, to calculate the loop diagrams contributions we will use the dimensional regularization method 
and the vanishing external legs approximation 0, Q which turns out to be legitimate because of the low masses and 
energy of the fermions with respect to the electro- weak (Mw) and SUSY (Msusy) breaking scales. When one does 
not use this approximation, the loop integrals to perform will have several different mass scales. Such calculations 
are much more complicated and need advanced mathematical methods as in [38|. In this framework, some specific 
functions will appear: 

- For the self-energies and the penguin: 



Ho(x,y) = ^fxy 



xhix 



_(x - y)(x - 1) 



+ {x <-> y) 



Go{x,y) 



(x-y)(x- 1) 



+ (x <-> y) 



For the box diagrams: 

H'(x,y,z) = y/xy 



xhix 



y]ny 



'AvlZ 



(x- y)(x- z)(x- 1) (y - x)(y - z)(y - 1) (z - x)(z - y)(z - 1) 



G'(x,y,z) 



y 2 lny 



(x - y)(x - z)(x - 1) (y - x)(y - z){y - 1) (z - x)(z - y)(z - 1) 



A. Vertices 

Following the notation of [H, we denote by JVy by \i = ^ij^j where ^ T = (B, W3, h u , ha, s), U and V are 
the 2x2 matrices required for the diagonalization of the chargino mass matrix. 



I't 



■ "V,rf^2 (1-75 
g V2M w cos,8 2 



m n4 v h (1+75 



FIG. 1: R-parity conserved vertices involving charginos and up fermions. 
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FIG. 2: R-parity conserved vertices involving charginos and down fermions. 
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FIG. 3: R-parity conserved vertices involving neutralinos, fermions and associated sfermions. 



We see in Fig.([T]) and ([2|) that vertices between chargino xt > a fermion / and the right scalar partner f R are 

m j 

My, 



negligible w.r.t. vertices between chargino \t ■> a fermion / and the left scalar partner fz because -j^- 1 for 



/ = /i, u, d. We will neglect the loops including the primary vertices below. 

The coupling of one neutralino x°l to one fermion and the associated left (right) scalar partner G J j L ^ is: 

G* fL = Q }S inO w N h + -A-N h 

where f = v,e, u, d. 

G fR = sign(m x o) 



QfsinOwNji + ° R n N j2 
cos9 w 



where / = e, it, d. 

The coupling of left (right) scalar partners to the Z°-boson is: 



C L(R) = T 3(/i(i?)) - Qfsin 2 b w . 
The coupling of scalar partners in their mass basis to the Z°-boson [4l[ is: 



i=l 

where tl = r^f^. The coupling of charginos to the Z°-boson [ill. |42|. is: 

0% = -VixV* x - \v i2 V* 2 + 6 ljS in 2 9 w 

and 

_ TT TT* }_TT TT* I X „„„2 

2 

We will assume below that all the parameters arc real 
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FIG. 4: R-parity conserved vertices with Z°-boson. 



B. Self-energy diagrams 

The computation of the supersymmetric contribution to Ap requires the evaluation of the Feynman diagrams. 
The self-energy, penguin and box contributions to Ap* can be written as 

We consider here the contributions to neutrino and antincutrino scattering involving self-energy diagrams. 
The first contribution implies a slepton-chargino loop (FigO^). We diagonalizc the stau mass matrix and use the 
vertices of Fig.fl]). We neglect the p,R contribution. 



where X a b 



V(s) = ZU u l pLi rf, [Go(x 4fk , 1) + 

Then, we have a neutralino-sneutrino loop (FigJBJ.). 



(13) 



l) + ln- 



- {v T 



(14) 



C. penguin diagrams 

We consider here the contributions to (anti)neutrino scattering involving penguins. We use the vertices of Fig.Q. 

Ap p (£) implies a chargino-slcpton loop in which the slcpton couples to the Z°-boson (Fig|5t). In Ap p (x + ), the 
chargino couples to the ^°-boson (Fig(5}y). Ap^{v) implies a neutralino-sneutrino loop in which the sneutrino couples 
to the Z°-boson (FigJBJ) and in Ap p (x°) t ne neutralino couples to Z°-boson (Fig[6p). 
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FIG. 5: R-parity conserved penguins and self-energies involving charginos and the charged Higgs. 
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FIG. 6: R-parity conserved penguins and self-energies involving neutralinos. 
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Ap£( x °) = frEl J= iGl L Gi L {N. l4 N j4 -N. l3 N 33 }x 

2H (X x o, tl , X x o 0tl ) + |g (X x o^ , X x o Ptl ) + ln^j - (P n -+ J 

We have neglected the loops involving ft,R. 

When we sum all the contributions (Fig[5]3, [5]5, [5£) involving the charged Higgs boson Q, we find: 
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FIG. 7: R-parity conserved box diagrams. 



The charge conjugation operators in the vertices of Fig.@ imply that all box diagrams involving charginos induce 
radiative corrections to neutrino scattering only, despite the fact that the neutrino fermionic lines of Fig.([7H) and 
Fif.j7]\) arc oriented on the left. As we did before, we diagonalize the stau mass matrix and we neglect the £ir 
contribution. The box diagrams involving charginos (resp. FigJTJc, [TJX, [7^) are: 



J2i=i r rL H '( X x^O^ X xt^ ,XfiP ' 

M OT (X+) = &PtoAX + ) fa - U L ) 



H '( X xt^ 



X xtv, 



vS] 



(19) 



(20) 



8?r mT 



- G '( X ^ L d L > & > x *£ & ) } • 



(21) 



The box diagrams involving ncutralinos with crossed fermionic lines will only contribute to neurino scattering 
because the neutrino fermionic lines are oriented on the right. Here, fn contributions are non-negligible. 
We consider here the box diagrams involving neutralinos for neutrino scattering (Fig[7J;): 



with left sfcrmions (e-L, ul, g?l) 
A/4(X°) = 



g w v^ 5 r< k ni c k ci j£a 

2tt l~,j,k=l '~ T vL'~ T vL' JT fL^ T fL 



fl. 



(22) 
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. with right sfcrmions (Sr, ur, cIr): 

/r l 



(23) 



Note that the contributions from all box diagrams involving antineutrino scattering are identical to the previ- 
ous contributions, only the forms of the boxes will be different. For instance if a ladder box was contributing to 
antineutrinos then the corresponding box for neutrinos will have crossed fermionic lines and vice-versa. 

E. Numerical results 

We make here a low-energy study where we take first generation sleptons degenerate with the second generation 
ones, and only allow the third generation sleptons to have a different mass. It is a possibility that sfermion masses 
may dynamically align along the directions, in flavour space, of the fermion masses, suppressing FCNC but allowing 
large mass splittings [43j . 
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FIG. 8: e as a function of p. We fix A = 0.4 and re = 0.5. Mi = 66GeV, M 2 = 133GeV and M 3 = 500GeV. The figure on the 
left represents the normal supersymmetric hierarchy (m.f = 300GeV and mp = 200GeV), the figure on the right represents the 
inverted supersymmetric hierarchy (m? = 200GeV and mp = 300GeV). 



We consider two experimentally allowed cases for the sleptons with a splitting between the third and the second 
generation: 

. mf = SOOGeV and = 2Q0GeV (normal supersymmetric hierarchy) 
. m,f = 200GeV and = SOOGeV (inverted supersymmetric hierarchy) 

We also assume that squarks are much heavier than sleptons: Mq = ITeV. These are the effects of gluino masses 
in the renormalization group evolution of scalar masses. 

The splitting among the sleptons of the second and third generations is mainly responsible for the size of Ap. 

The purpose of this Section is the investigation of the supersymmetric parameter space in some specific cases. To 
this end we made a subroutine to the Fortran code NMHDECAY, which is available on the NMSSMTools web page 
j4p] |. |44j |. [4o| . This subroutine computes the different R-parity conserved supersymmetric contributions to e. Note 
that e is the same for antineutrinos, as in the SM. 

By making scans over the supersymmetric parameter space, we can obtain, in some regions of the parameter space, 
divergences for either m x + = rrif k , m x + = m p, L or m x° = m ^t • This is because we assume vanishing external legs. 

In supersymmetry, there are many parameters, e doesn't depend very much on A and re so we fix them as we usually 
do in supersymmetry: A = 0.4 and re = 0.5 . We allow some other parameters to vary: tan/3, Ma and /i e ff- 

tan/3 is the ratio v u /vd where v u and Vd are the v.e.v. of the scalar Higgs fields H u and Hd- Ma is an effective 
supersymmetric parameter somewhat equivalent to the second pseudoscalar Higgs mass in our regions of the 
parameter space. 
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FIG. 9: e as a function of /i. We fix A = 0.4 and k = 0.5. Mi = 150GeV, M 2 = 300GeV and M 3 = ITeV. The figure on the 
left represents the normal supersymmetric hierarchy, the figure on the right represents the inverted supersymmetric hierarchy. 



We show two different illustrative situations motivated by high-energy models for the gaugino masses: 



Mi = 66GeV, M 2 = 133GeV and M 3 = 500GeV (cf Fig.©) 



Mi 



lbOGeV, M 2 = SOOGeV and M 3 = ITeV (cf Fig.®) 



In Fig.©, tan/3 varies between 2 and 15, M^ varies between 579 and 2000 GeV, \i varies between 300 and 800 
GeV. We see in Fig.© that contrary to the SM case, V^ r can be either positive or negative. 

In Fig.©, tan/3 varies between 2 and 18 (for the figure on the left) and between 2 and 9.6 (for the figure on the 
right), Ma varies between 500 and 1000 GeV, \i varies between 200 and 564 GeV. 

On the other hand, as we can see in the Fig.©, e can go up to 2 x 10~ 2 in this region of the parameter space. The 
sign of \i does not have an important impact on the maximal value that e can reach. 

We finally mention that other extensions of the SM may also lead to sizable effects upon An Tfl . In the next section, 
we will consider supersymmetric i?-parity violating interactions which can have important effects on the neutrino 
indices of refraction already at the tree-level [46[ . 

IV. R-PARITY BREAKING SUSY CORRECTIONS 
A. Introduction 

The superpotential Wnmssm is not the most general supcrpotential we wan write because there exist some other 
gauge invariant couplings that we didn't take into account. These new couplings break a discrete symmetry called 
R-parity. This symmetry requires that an interaction must have an even number of SUSY particles. If R-parity is 
broken, the LSP will no longer be stable [4?| . R-parity violating interactions also violate lepton number (L) or baryon 
number (B). To be as general as possible, the superpotential has to contain the following terms: 



W fip = (~2^ykLiLjE^ + \' ijk LiQjDy. + —X'^apj U^D^D^ 1 + Li i H u L l + XiSH u L^j 



The lagrangian Cfc p can be derived using the common procedure |37l . |48 
If we consider this part of the superpotential: 



masses 



Wnmssm + fJ>iH u Li + XiSH u Li 



(24) 



(25) 



we can generate two neutrino masses at tree- level [37J by mixing neutrinos and neutralinos. By considering a 
See-Saw-like mechanism, the mass matrix is: 
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M- = [ M ™ SSM £L ) (26) 



^ p 3X 3 



where Mnmssm is the R-parity conserved neutralino mass matrix and £^ is the part of the mass matrix induced 
by R-parity violation which mix neutralinos and neutrinos. We will assume here: Wj/wud *C 1, l/^i/VI <C 1 and 

<C 1 in order to reproduce the neutrino phenomenology. Vi are the v.e.v. of the sneutrinos. 
This model is self-consistent because we give here a way to generate the neutrino masses contrary to many other 
scenarios of radiative corrections on neutrino indices of refraction. 

In the following, i/j = Nij^j where ^> T = (B, Wa, h u , hd, s, v e , v^, v T ), i = 1..3 and j = 1..8. v e , and v T represent 
the neutrinos in their flavour basis, Vi represent the neutrinos in their mass basis. Xi = ^(i+3)j^j where i = 1..5 and 
j = 1..8. x° are the five NMSSM's neutralinos. 



B. Tree-level corrections 



Because of R-parity breaking, new interactions are possible between SUSY particles. This yields new one-loop 
corrections but also tree level corrections contrary to the R-parity conserved case. Such tree level corrections in 
the Supernova context have been partially studied in prl l49j . Since this work focuses on V^ T , we only concen- 
trate on graphs involving fi and r neutrinos. Moreover, we only consider neutrinos which do not change of flavour 
after the R-parity breaking interaction, therefore we are not concerned by off-diagonal term in the matter Hamiltonian. 




FIG. 10: R-parity broken Tree-level corrections. 



Here, we give a specific contribution to antincutrino scattering (TigflOk) with left slcptons i^.i- 

A Pa (l L . z ) = - V (MiLZ^k) M w (27) 

With left down squarks dtA (FigflOb): 

&Pb(d L ,i) = Ap a (£ Lti ) (Afeii -* X' kil ,£L,i -> d L ,i) (28) 

The right slcptons Ira and the right down squarks cLr^ do not contribute. 

The following contributions induce corrections to both neutrino and antincutrino scattering: 
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for diagram (FigfTUb): 

5 

Ap c (x°) = - (%+3)8^(i+3)8 ~ N u+3)7 N {l+3)7 ) {N u+3)4 N {l+3)4 - N u+3)3 N {l+3)3 ) (29) 
3,k=l 

for diagram (FigfTUh): 

Ap d (4,i) = Ap a (i L ,i) (30) 

for diagram (FigfTOe): 

Ap e (d Lii ) = Ap b (d Lli ) (31) 
For (FigfTUE) and (FigflOb). the right slcptons and the right down squarks dn,i do not contribute. 

C. Self-energy corrections 

We willingly leave certain factors not simplified to be able to compare with other corrections more rapidly. All the 
contributions will be equal for neutrino and antineutrino scattering. 





diA i 
s z 
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FIG. 11: R-parity broken Self-energy corrections. 



For Fig.CHJ): 



3 (\2 _ \2 \ ( 771 ~ ~1 

V* (E) = £ \ Go(X, Jrj , 1) + ln-^ I , (32) 



a/"(e) = a/«(e) (A Wi -> a;.,,.'/,.., - 4,,-) 



For Fig.((nh): 



Ap iL (E) = A/*(E) (A Wi -» A^,^ - l Lj ), (33) 
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In the R-parity breaking scenario, there are NMSSM specific self-energies due to the term \iSH u Li in the superpo- 
tcntial. This term induce one new coupling between a neutrino, one Higgs (scalar or pseudoscalar) and a neutralino. 
We present here the contributions for scalar Higgses and for pseudoscalar Higgses. 

For both Fig. (fTTk ) and (fill ), the contribution from scalar higgses: 

V(S) = EB« + ^ + 3)3^) [Go(X x o hj , 1) + m^j (34) 

from pseudoscalar higgses: 

AP Q (S) = + ^ +3)3 ^) [Go(X x o aj , 1) + In^j (35) 



D. penguin type diagrams 



( L,RJ 




d L,R,j 



di 

z f 



V 


k V 










/ 


<, z f 




di v 










f 


X z f 



FIG. 12: R-parity broken penguin type diagrams. 



We consider here the penguin diagrams. All the contributions will be equal for neutrino and antineutrino scattering. 
For the diagram of Fig.fHJ) and Fig.(M): 



Ap e -(£) 



aw 
16tt 



3 (A^-A 2 



E 



2 jit 



G (X ir l)+ln(^i) 



(36) 



Ap^(d) = Ap^(£) (X kji -> X' kji ,I Ltj -> dij^ 
For the diagram of Fig. (fT3k ) and Fig.(fl3hi): 

Vp(^) = - 7 -c1 ^ 



G (X^.,l)+ln(™^) 



(37) 



ApJ(di) = Ap e P (t L ) (Xkji -> A' fej4 ,^j -> di,,-,ii -> di) 



For all the previous penguin diagrams, the right sleptons £rj and the right squarks dn,j do not contribute. 

We have below NMSSM specific penguins also due to the term \iSH u Li in the superpotential. 
For the diagram of Fig. (TT3"b ): 



A/£(x°) 



8tt 



^(j+3)4A r ( l+3 )4 



^(j+3)3^(i+3)3 



2i?o(^ x o/ lm , A x o, lm ) + (G (X x o hm , X x o hm ) + ln^p-} 



(38) 



15 



v On v v Xi v 




FIG. 13: R-parity broken NMSSM specific penguin type diagrams. 
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FIG. 14: R-parity broken box diagrams. 



For the diagram of Fig. ([l~3"h ): 



a (v°) — q u , A 3 -A 2 



zCi.j = l zCn=l -^(l 



,r r ^1,-7=1 ^n=l- (i+3)5-^(j+3)5-Pnl (%+3)4^(i+3)4 - N( j+S)3 N( i+3)S ) 



We have two identical contributions from the diagram of Fig.(fl~3"b) and Fig.(fl~5b) and we obtain: 



aw A? — A 



2 5 3 2 
3 /x 2 



4tt g 2 



^2 X! ^+3)5^ml P ™l^M™ ^ G o( X h mX °' X a n X°) + U1 ~ 



i— 1 m— 1 ri—1 

where = SVni-Pni — £m2-Pn2 following the notation of [io| . 



(39) 



(40) 



E. Box diagrams 

i. R-parity broken box diagrams 



We present in this section the box diagrams corrections. Each box diagram correction involves Fierz transformations 
to obtain a similar form as the SM tree level interaction. The box diagrams with crossed fermions lines (Fig |14b 2. 
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FigfT4l2. FigfT4fe2. FigfT4h2. Figfnj^, FigU4W2) will only contribute to ncurinos, meanwhile the uncrossed fcrmion 
lines (TigfTlbl, FigdU-l, Fig[T3fcl, FigfTHil, FigfTiVl. Fig[T4"Wl) will contribute to antineutrinos. However, the 
contributions to antineutrino scattering will be exactly the same as the contributions to neutrino scattering. We show 
below the contributions to neutrino scattering. 
For the diagram of Fig. ([i"4V 2) with t^^ — d R:Tn : 

A v2 in t \ a w Mw ( A 3ij X 3i k - A 2u - A 2 »fc)A^ lm A / fclm 

Ap box (£ L -d R ) = - — 2^ ^2 Ja G \ x t j M^ X i*,d R , m i x i L , i ,dR,J- ( 41 ) 

i,j,k.m—l dfi :m 

For the diagram of Fig. (JT3J:2) with d,L,i — (-L.m an d the diagram of Fig. (JT5J>2) with i&n — «l m : 

Ap r b 2 ox (d L - £l) = Ap v box (£ L - d R ) (A -> A', \' jlm \' klm -> A^A^,^ -> dj,k,d R , m -> ^l,™), (42) 



Ap s b 2 ox {d L - 2i) - Ap r b 2 ox {d L - £ L ) {X' mlj X' mlk -> Ai mj Ai mA .,£i, m -> Si, m ). 
For the diagram of Fig. (jT3t2) with — ^L,m and the diagram of Fig. (Ti~4"k 2) with — ul. to : 

i,j.k,m— 1 &L,m 

Ap S b l x (d R - U L ) = Apll x (d R - £ L ) {KnljKnlk -> XimjKmk>^L,m ~* "i,m)- 

For the diagram of Fig. p^b 2) with — VL.m and for the diagram of Fig. ([i"4"u 2') with £ Ri i — i>L,m' 

^f ox (d R -H) = ^ f (^A - ^^^^'^'//-.V,,, ...V,,.,,,...*,,,, „, ). (44) 

i,j,k,m=l v L,m " 

A Pbox(^R ~ vl) = Apf ox (d R - v L ) (KnjiKnki x 'mij A ™i fe thenA' -> A,d,-, fc -> £j,k,d R ,i -> £ R ,i). 
For the diagram of Fig. ([1^ 2) with — d R m and for the diagram of Fig.(fT4"W2) with li,^ — UL, m '- 

Ap box {l R -d R )-— ^ ^2 -A ^^^'^^'^V)' ( 45 ) 

i,j,k,m=l dn jJrl 

Apl"o X {£L - U L ) = Ap% ox (l R - d R ) {X' jlm \' klm — > X' jrnl \' kml , d R>m — > UL, m ,tR,i ~ * £l,i)- 

For the diagram of Fig. (fl4T 2) with c?L,i — d R>m : 

A^2 /J 7 n _ a W \^ M W ( X 3ij X 3ik - X' 2 ijX' 2ik )\" mj \'{ mk „, , y y y x 

^Pbox, a *{dL-d R ) - — 2^ 2^ ^2 72 e<T 7 5 e "7/3 H ( A d|,< m ' A 4,dl m > X d L ,i,dlJ- 

i,j,k,m=l ,8,5,7=1 rf^ m ^ 

(46) 

For the diagram of Fig.|l4^2) with <iR : ; — u R , m : 

3 3 , r 2 ('X' \' _ \' \' \\" \" 

^Pbox,0<r\"'R~ u R) — 2^ 2^ rrr; 72 e aSa e ay0 H (A d j )figm , a^^, a^^). 

i,j,fe,m=l a, 5,7=1 u R,m 

(47) 
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The indices of e denote SU(3) color indices where e is antisymmetric, e Q /3 7 = = e 7Q ^. For example, for the 
diagram of Fig.(|14b2) with dn — Uu >m , the down quark on the left carries the index (3 and the down quark on the 
right carries the index a: 



M^r (^3ij^3ik ^2ij^<2ik)^mlj^mlk 



(48) 



i,j,k : m= 1 a, (5,7— 1 u h 

H '( X d°U 



The diagram of Fig. (|14b 2) with G^.j — m can be deducted from the previous one by removing the e factors and 
by taking a factor 4 and A" — > A', — > z>l: 



% 

47T 



E 



(A3j,A 3ife A 2 j,A 2ifc )A m ,iA TOi:1 



i.j,k,m— 1 



(49) 



Finally, for the diagram of Fig. (|14h 2) with , — VL,m'- 

^pfoAh - H) = Apf ox (d L - v L ) (A' A, dj, k -> ^, fc , d il4 -> fx.i). (50) 
The rest of diagrams do not contribute because of a different form as the correct tree level SM form given in Eq.([2]). 



2. Cancelling diagrams 
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u / ^mi e m N u 




V 


k 
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v H v 
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di 


V 
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h,k,'\hj 

d , ' u m - N d 


z2) 


d-R,k\ 
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e / u m \ e 
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FIG. 15: R-parity broken radiative corrections which cancel out. The Feynman graphs with uncrossed scalar lines on the left 
cancel with the Feynman graphs with crossed scalar lines on the right. 



It is interesting to notice that amongst the possible radiative corrections induced by the R-parity breaking inter- 
actions, some of them cancel. Contrary to the previous box diagrams, all box diagrams in Fig. (|15[) should contribute 
equally to neutrinos and anti-neutrinos because the fcrmionic lines do not link the neutrinos to the fermion present 
in the external leg. The particles exchanged between the two fermionic lines here are all scalar. 
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To understand the consequences, we write the tensorial part of the scattering amplitude of a ladder scalar line box 
diagram (FigfTSkl) and the corresponding crossed scalar line box diagram (Figll5k2): 
For the left diagram we have: 



Mi = u(ki) 



dq 4 .. 1 + 75^ + ^) ., 1-75 , s 
»A 3 fei — 2 —i-^ji — 7. — u (Pi) 



u(k 2 )iX mk _ 
while for the right diagram we have 

M 2 = u{ki) 



(27T) 

1 - 75 i($ + m Vm ) l+ 75 

z A 



(pi - q) 2 - ™% . - i) 2 ~ m i 



q 2 — m 2 , 



mjl- 



u{p 2 ) 



dq 4 .. l + 75^ + m O.. 1-75 , s 



(27T) 



(pi — q) 2 — m~ (fci — q) 2 — m 



u(k 2 )i\ m ki 



1-75 itife - Vi + i + ™>u m ) . . 1 + 75 



9 9 

q- - mi 



i A 



"(Pa) 



(51) 



(52) 



Therefore , in the approximation of zero external legs (pi, p 2 , ki, k 2 = 0), we have Mi = — M 2 . From a physical 
point of view, it is actually quite natural that such corrections cancel out when considering vanishing external legs, 
the Fcynman box diagrams are anti-symmetric in this case since scalar particles are exchanged and the sign of the v m 
neutrino impulsion is opposite. 



V. IMPLICATIONS ON THE SUPERNOVA NEUTRINO FLUXES 



In the supernova environment the density is sufficiently high that neutrinos encounter not only the so called high 
resonance associated with #13 and the low resonance associated with Bi 2 but also the \i — r resonance whose precise 
conditions depend on the hierarchy 50] . The importance of such correction term has been investigated in (2{| [2?| 
showing that such resonance could influence the electron neutrino flux. This radiative correction term may also 
influence the electron (anti-) neutrino in an indirect way. Indeed, V^ T breaks the symmetry between muon neutrinos 
and tau neutrinos and creates a CP-violation dependence for the electron (anti-) neutrino surviv al p robability and 
consequently on the electron (anti-) neutrino flux. Such a phenomenon was numerically observed in [5 II ] - The addition 
of the non-linear neutrino-neutrino interaction induces larger effects on the electron neutrino fluxes up to a level of 
10% inside the supernova (29[. In this section, we analytically demonstrate that the inclusion of V^ T (or more generally 
the inclusion of an interaction term in the total Hamiltonian breaking the symmetry between mu and tau neutrinos) 
will make the electron survival probability dependent upon the CP-violation phase S. The fact that in the SUSY 
framework, V^ r can be up to a factor 2 x 10~ 2 V e implies sizeable effects on the electron neutrino fluxes seen in a 
detector on Earth. 



1. The factorization 



In a dense environment the neutrino evolution equations with one-loop correction to matter interactions are given 
by: 

( % ) 

where ^ a denotes a neutrino in a flavour state a =e, fi, r, £^=1.2, 3 being the energies of the neutrino mass eigenstates, 
and U the unitary Maki-Nakagawa-Sakata-Pontecorvo matrix 

(I \ / ci 3 si 3 e~ iS \ ( c 12 s 12 0\ 
U = T 23 T 13 T 12 = c 23 S23 1 -S12 c 12 , (54) 

V -S23 c 23 J \ -sis e 18 da / \ 1/ 

where c.y = cosOij (s.y = sin9jj) with #i 2 ,# 23 and 613 the three neutrino mixing angles. The presence of a Dirac 
6 phase in Eq. (|54p renders U complex and introduces a difference between matter and anti-matter. The neutrino 




U 



Si 
E 2 
E 3 




19 



interaction with matter is taken into account through an effective Hamiltonian which corresponds, at tree level, to the 
diagonal matrix H m = diag(V c , 0, 0), where the V c (x) = y/2GFN e (x) potential, due to the charged-current interaction, 
depends upon the electron density N e (x) (note that the neutral current interaction introduces an overall phase only). 
Following the derivation of [1^, [5l[ , to obtain explicit relations between probabilities and the CP- violating phase, it 
is convenient to work within a new basis where a rotation by T23 is performed. In this basis, one can factorize the S 
matrix, defined by diag(l, 1, e ), out of the Hamiltonian, so that: 





rpf rpO t 

± 12 J 13 



(55) 





C23S23e 4<5 Vw 



,2 y 

'(IT 

-C23S23e "° VfiT 



Con Vn 




Contrary to the case where only tree level matter interaction is considered, the S matrix does not commute with the 
matter Hamiltonian in the T23 basis. This fact implies that: 



and, therefore, that 



H T (S)^SH T (S = 0)S\ 
U m (8) + SU m (8 = 0)5t. 



(56) 
(57) 



2. Consequence on the electron neutrino survival probability 

Nevertheless, the factorization of the S matrices is always possible but the Hamiltonian H't{8) will depend on 8. 
We can rewrite Eq.Q55|) in the evolution operator formalism to be: 

A eh A m A^ 5 = H' T {8) A eh A m A? fi e iS | (58) 
A ~ P -' iS A~~p~ iS A-~ I \ A ~ P ~ iS A~~p~ iS A~~ 




Defining the Hamiltonian H't{8) by 



2. 



b c 
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H' T {8)= ( b d (e-ge ld ) | (59) 

c (e-ge-* s ) f 

we can rewrite the evolution equations for the amplitudes of the first column of the evolution operator, which corre- 
sponds to the creation of an electron neutrino v e initially: 

i——A ee — a A ee + b A e n + c A e fe 
dt 



i-rAefi = bA ee + dA ejl + (e-ge z5 )A e fe tS 



d 
~dt 

ij t Aefe- tS - cA ee + (e-ge u )A eFl +fA ef e- lS (60) 

Similarly, we can write the same equation for the amplitudes B a f3 when 8 is taken to be zero and look at the difference 
between the amplitudes that depend on 8 and those which do not. The basic idea here is to prove that, because of the 



2 Note that the terms a, b, c, d, e ,f and g are real. 
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one- loop correction term V^ r , the function (A ee — B ee ) can not remain zero function by showing that its derivative is 
non zero. 

(61) 



l Jt {Aee 


— B ee ) 


= a (A ee 


— B ee ) ~ 


f b (A ejl 


- Befi) " 




— B e f) 




.d 


- B ep) 


= b(A ee 


- B ee ) ~ 


r d {A efl 


~ Bey) " 


V e (A ef e- lS 


~ B e f) " 


Vg{A e f-Bef) 




— B e f) 


= c (A ee 


— B ee ) ~ 


h e (Aef, 


~ ^e A ) " 


h / (A ef e-* s 


— B e f) - 


- 9 (A e jie~ i5 — B e jj) 



d 
l dt 

Let us now take a closer look at Eqs. (|6ip . The initial condition we are interested in, namely a v e created initially 
means that initially the amplitudes A and B are: 

| = | | =H A efl I = I B efl I = I | (62) 

When <7 = (V^ r = 0), it is easy to see that inserting the initial conditions into Eqs. lfBTj) will imply that the functions 
f e = A ee — B ee , ffj. = A e p, — B e p\ and f T = A e fe~ lS — B e f will be equal to zero. By discretizing time we see, by 
recurrence, that if those functions are zero at beginning, they will be equal to zero at all time. But when g ^ 
(V^t 0) we have to look also at the evolution of the functions = A e jj,e^ lS — Be, and f T = A e f — B e f- Their 
respective evolution equation can be easily derived from Eq. (|60p to yield: 




A e e 




'Bee 


A e jl 


H 


B e [i 


Aefe~' lS 




y B e f 



4(^e-« - Be,) = b (A ee e-- - B ee ) + d (A^ - B*) + e (A^™ - B e ,) + , (^e"« - B ef ) (63) 



ij t (Aef - Bef) = c (A ee e lS - B ee ) + e (A e ,e lS - A e jx) + f (A ef -B ef )-g {A e[l - B e ,) (64) 



i-(4M*e-" -B e ,)(t = 0) = i-U(t = 0) = b(e 



Initially, the derivatives are : 

-id 

■J[i{t. = <J) = uyc 

lj t (A e .f - Bef)(t = 0) = ijJ T (t = 0) = c(e lS - 1) (65) 

We just proved that since the functions / M and f T are non constant zero functions, the functions / M and f T won't be 
zero as well. But does it implies that the function f e is non zero at all time? No, because the contributions from / M 
and f T could cancel in the evolution equation (|6~lj) of f e . To precisely study the evolution of / e , we discretize time 
such as t = N * At with N G N and 

d feit + At)-Ut) 

dt J At y ' 



Using Eqs. lfBTj) and the time discretization we see that: At t = At: 

i 

c(e iS - 1) 
MM) = { . ' At 



which implies that: At t = 2At: 

iS i\a2j 



(67) 



/„(2At) = gb{e- M -1)AH 
f T (2 At) = -gc(e l& - l)A 2 t 

(68) 
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leading at t = 3At 

/ e (3Ai) = t (<?6c(e- t5 - 1) - gbc(e lS - 1)) A 3 t 

= ^gbc(e~ l5 - e lS )A 3 t 

= -2gbcsm5A 3 t (69) 

This last formula proves that the function / e is not the constant zero function when 6 ^ 3 , therefore A ee ^ B ee and 
consequently, for V^ T ^ 0: 

p(i/ e -► i/ e , s + o) ^ p(^ e -» i/ e , <y = o) (70) 

Therefore, when 5 is non zero, it has an influence on the value of P(y e — > v ei S). The luminosity of a neutrino emitted 
initially as a flavour a is 

L v Jr, E v ) - T sj Ev j F2{ri) 1 + exp {Ev j Tvjt _ v) 

where F2{rj) is the Fermi integral, L° Ua and are the luminosity and temperature at the neutrinosphere. The v e 
and v e fluxes will depend on 8 even when the luminosities L v and L Vt are taken equal at the neutrino-sphere: 

<j> Ve (5) = L^P(v e -^v e ,5) + L VjL -► i/ e ) + P(i/ T v e )) 

= L VsL P(v e -> i^ e , (5) + (1 - P(v e -> i/ e , (5)) 

= (i^-^jP^^^^ + L^ (72) 

This analytical derivation proves that, if the dependence on 5 of the evolution operator cannot be factorizcd then the 
electron neutrino survival probability depends on 5. Nevertheless, this implication had been observed numerically. 
We can easily generalize this derivation to other interactions that would distinguish from v T like for instance, 
non-standard neutrino interaction (52j. Therefore we can state that as soon as the medium effect on and on v r is 
not the same, effects of the CP-violating phase on the electron neutrino (and anti-neutrino) fluxes will appear. 



VI. CONCLUSIONS 



In this paper, we have investigated the radiative correction on the fi — T neutrino indices of refraction coming from 
beyond standard physics. In the NMSSM, we have shown that the sign of V^ T depends upon the hierarchy of the 
sleptons masses and therefore could be negative contrary to the Standard Model case. After writing and adding a 
subroutine to a low-energy code taking into account all current constraints on SUSY we showed that e can increase 
up to the order of 2 x 10 -2 depending on the supersymmetric parameters. Such value could be highly important in 
the calculation of neutrino fluxes from core-collapse supernovae as it can induce sizeable effects upon the electron 
(anti)-ncutrino fluxes. In a second part we have calculated all contributions from It-parity breaking interactions on 
the radiative corrections V^ T . Taking into account such interactions, we showed that NMSSM distinguishes from 
MSSM in this case and bring new possible contributions. The next step in this type of calculations would be to 
see the consequences of theses values of V^ T on the supernova neutrino fluxes and to use it in order to survey the 
supersymmetric parameter space (53|- Secondly, we would have to calculate corrections with gravitino loops and all 
contributions from R-parity breaking interactions. Another interesting possibility would be to calculate the radiative 
corrections for the neutrino-neutrino interaction. Such calculation has been recently done in the SM [54[ and SUSY 
framework could yield potentially much bigger effect. Finally, as an application for V^ lT wc demonstrated that the 
inclusion of such a term implies that the electron (anti-) neutrino survival probability, and consequently the electron 
(anti-) neutrino fluxes, will depend upon 8. The consequences of such dependence and the fact that V^ lT can be up to 
2 x 10~ 2 V e will be studied in a future work. 



3 Note that with such formula, f e is also equal to zero when <5 = it, but going to the forth step will show, after a tedious but straightforward 
calculation, that f e is non zero even when 5 = it. 
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